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On the progenitors of millisecond pulsars by the recycling 
evolutionary channel 
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ABSTRACT 

The recycling model suggested that low-mass X-ray binaries (LMXBs) could evolve 
into binary millisecond pulsars (BMSPs). In this work, we attempt to investigate the 
progenitor properties of BMSPs formed by the recycling evolutionary channel, and 
if sub-millisecond pulsars can be produced by this channel. Using Eggleton's stellar 
evolution code, considering that the dead pulsars can be spun up to a short spin 
period by the accreting material and angular momentum from the donor star, we 
have calculated the evolution of close binaries consisting of a neutron star and a 
low-mass main-sequence donor star, and the spin evolution of NSs. In calculation, 
some physical process such as the thermal and viscous instability of a accretion disk, 
propeller effect, and magnetic braking are included. Our calculated results indicate 
that, all LMXBs with a low-mass donor star of 1.0 - 2.0 M Q and a short orbital 
period (< 3 — 4d) can form millisecond pulsars with a spin period less than 10 ms. 
However, it is difficult to produce sub-millisecond pulsars by this evolutionary channel. 
In addition, our evolutionary scenario cannot account for the existence of BMSPs with 
a long orbital period (P rb > 70 — 80d). 
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1 INTRODUCTION 

Millisecond pulsars (MSPs) and normal pulsars have dis- 
tinct observed properties, and they inhabit two different 
regions in magnetic fiel d - spin period (B — P) diagram 
l|Manchester et al] [2005). Normal pulsars have a spin pe- 
riod of P ~ 1 s and a magnetic field of B ~ 10 12 G. How- 
ever, MSPs show some distinct observed properties such 
as short spin period (P < 20 ms), low spin-down rate 
(P ~ 10~ 19 — 10~ 21 ss" 1 ), old characteristic age (r = 
Pj (2P) ~ 10 9 — 10 10 yr), and weak surface magne tic fields 
(B ~ 10 s - 10 9 G) (|Manchestedl2004l : lLorimerll200SD . About 
75% MSPs are in binary system (called binary millisecond 
pulsars, BMSPs), whereas that is only < 1% for normal pul- 
sars. 

At present, there exist two scenarios to account for the 
formation of MSPs. The first one is the recycling model, 
in which MSPs are proposed to be the evolutionary prod- 
uct of neutron star (NS) low-mass X-ray binari es (LMXBs) 
or in termediate-mass X-ray binaries (IMXBs) l|Alpar et al.l 
1982). The pulsar crossed the so-called deathline accretes 
the mass and angular momentum from the donor star that 



overflows its Roche lobe, a nd can be subsequently spun up to 
a mi l lisecond spin-period ( Bhattachar va fc van den Heuve] 
Il99ll :l Tauris fc van den Heuvelll2006l ). During accretion, the 
magnetic field of the NS decrease to be B ~ 10 8 — 10 9 G due 
to a ccretion-induced field decay (|Konar fc Bhattacharval 
Il997h . When the mass transfer ceases, a BMSP consisting 
of a recycling NS and a low-mass (< 0.4Mq) helium white 
dwarf is produced. The discovery of the accreting millisecond 
X-ray pulsar Sax J 1808.4-3658 presented str ong support to 
this scenario |Wiinands fc van der Klislll998l ). Recent opti- 
cal observations also confirm that there exists a transition 
link between X-ray pu lsar and millisecond radio pulsar (e.g. 
I Archibald et al.ll2009r i. 

In another evolutionary channel, MSPs may be formed 
by accretion-i nduced collapse (AIC) of ONeMg white dwarfs 
l|Michel|[l987h . When the mass of an ONeMg white dwarf 
reaches the Chandrasekhar mass limit by accreting from its 
donor star, the electron-capture process leads to a gravita- 
tional collapse rather th an a Type la explosion, and results 
in the formation an NS (|Nomoto fc Konddll991l) Q. If MSPs 
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1 When the ONeMg core of an asymptotic gi ant branch star 
(ISiessll2007l; iPoelarends et alj|2008h or a He star iNomoto ||l987ft 
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formed by t he collapse of low fi eld (10 3 — 10 4 G) white dwarf 
population (| Jordan et al.ll2007h . their magnetic field should 
be in the range 10 s — 10 10 G, and without invoking sig- 
nificant field decay. Recently, the calculated birthrates by 
population synthesis approach indicate that the AIC chan- 
nel may play an important role in forming MSPs (Hurley et 
al. 2010). As an alternative formation of MSPs, this evolu- 
tionary channel has been widely explored by some authors 
ilWickramasinghe et aHl2009l;lHurlev et al.ll2010l ; IChen et al.1 



l2011al ). In particular. |Pu et all |2009l) argued that AIC pro- 
cess of massive white dwarfs can produce sub-millisecond 
pulsars (quark stars) with a spin-period less than 1 ms (or 
less than 0.5 ms). 

The purpose of this paper is to systemically explore the 
initial parameter space of LMXBs that could evolve into 
BMSPs via the recycling evolutionary channel. In addition, 
we also attempt to examine if this channel can form the so- 
called sub-millisecond pulsar. The structure of this paper is 
as follows. We describe the input physics that is necessary 
in the evolution calculation of LMXBs in section 2. The 
calculated results are presented in section 3. Finally, we give 
a brief discussion and summary in section 4. 



IChen fc Panel i|201ll) found that accretion disk instability 
model successfully reproduces the orbital period and the 
mass of the WD of PSR J1713+0747. 

When the mass tran sfer rate — Md is lower than the 
critical mass-transfer rate (Ivan Paradiislll996l ; lDubus et al. I 
1999) 



M c; 
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where P or b is the orbital period of the binary, the NS accretes 
only during outbursts. Defining a duty cycle d to be the 
ratio of the outburst timescale to the recurrence time0, the 
accretion rate of the NS Af ac = —Md/d. Otherwise for a high 
mass transfer rate —Md > M cr , we assume M ac = — Md- 
Certainly, the mass growth rate of the NS should suffer the 
limitation of the Eddington accretion rate (MEdd ~ 1-5 x 
1O _8 M0). The excess material is assumed to be expelled 
from the vicinity of the NS by radiation pressure, and carries 
away the specific orbital angular momentum of the NS. 



2 INPUT PHYSICS 

Using a stellar evolution code developed by Eggleton 
(see Eggleton 1971, 1972, 1973), which has been updated 
with the latest input physics over the p ast three decades 
|Han et al.Hl994l ; iPols et al. Ill995l . I199ST ). we calculate the 
evolution of binaries consisting of a NS (of mass Mns) and 
a main-sequence donor star (of mass Md)0, and test if they 
can e volve into MSPs. The stel lar O PAL opacities was taken 
from Rogers fc Iglesias I l| 19921 ) and I Alexander fe Ferguson] 
(|l994 ) for a low temperature. In our calculation, the ratio 
of mixing length to local pressure scale height was set to 
be 2.0, and the overshooting parameter of the donor star 
(with a solar chemical com position X = 0.70 , Y = 0.28, and 
Z = 0.02) is taken to be (|Dewi et al.ll2002i ). 



2.1 Accretion disk instability 

With nuclear evolution, the donor star overflows its Roche 
lobe, and transfer hydrogen-rich material onto the NS. Due 
to the high angular momentum, the accreting material forms 
a disk surrounding the NS. If the effective temperature in 
the accretion disk is below ~ 6500 K (the hydrogen ion- 
ization temperature), the disk accretion shou l d be ther- 
mally and viscous unstable l|van Paradiis l ll996l ;l King et al.l 
1 19971 ; lLasotal 120011 ). Meanwhile, the accreting NS will be a 
transient X-ray source, which appears as short-lived out- 
bursts phase and long-term quiescence phase. Recently, 



grow to a cr itical mass, electron-captur e supernova can also pro- 
duce an NS. iPods iadlowski et al. ( 2004) suggested that the min- 
imum mass of the NS progenitor may be 10 — 12 Mq for single 
stars, while this value can be 6 — 8 Mq in binaries. 
2 Certainly, NSs may also recycled by accreting the material 
from the He star companion. However, some studies show that 
the evolution products of NS + He star systems should be 
intermediate-mass binary pulsars or high-m ass binary pulsars (see 
iFrancischelli et aljliool IChen et alj|2011bl ) . 



2.2 Magnetic braking 

Low-mass donor star would be braked to spin down by 
the co upling between the ma gnetic field and the stellar 
winds (|Verbun fc Zwaan l ll98ll ). However, the tidal interac- 
tion between the donor star and the NS would continuously 
spin the star bac k up co-rotation with the orbital rotation 
l|Patterson Ill984l ). Therefore, magnetic braking mechanism 
indirectly carries away the orbital angular momentum of bi- 
naries. 

For the angular m o mentu m loss rate via magnetic brak- 
ing, iRappaport et al.l (1 19831 ) developed an empirical for- 
mula, i. e. 



-3.8 x 10 



■ 3O M 2j R|( J R 2 /i?0) 7 a; 3 dyncm, (2) 



where R2 is the radius, uu the angular velocity of the donor 
star, and 7 is a dimensionless parameter in the range of zero 
to four. This standard magnetic braking model is widely 
applied in studying the evolution of cataclysmic variables. 
However, studies on rapidly rotating low-mass stars with a 
spin period below 2.5 - 5 days in young open clusters show 
that the stand ard model overestimates the angular mom en- 
tum loss rate l|Queloz et al. I ll998l; lAndronov et~aT1l2003r ). 

In calculation, we adopt an indu ced magnetic braking 
description given bv lSills et all (2000), in which the angular 
momentum loss rate is 



J mb 



-Kuj d 



Rd M e \ 
R e M d J 



1/2 



1/2 



OJ > CJ cr it 



(3) 



where K = 2.7 x 10 gem |Andronov et al.ll2003h . Wcrit IS 
the critical angular velocity at which the angular momentum 
loss rate reaches a saturated state, ui = 2n/P OI ^ and Rd 
are the ang ular velocity and the ra dius of the donor star, 
respectively. iKim fc Demarauel (|l996l ) proposed that cj cr it is 



3 iKing et al] l|2003t ) proposed that the typical value of duty cycle 
is about 0.1 to 0.001. In this work, we take d = 0.01. 
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inversely proportional to the convective turbulent timescale 
of the star when its age is 200 Myr , i. e. 



r 

Merit — Wcrit Q , 

T 



(4) 



where £j cr it,0 = 2.9 x 10 -5 Hz, r©, and r are the convective 
turbulent timescales of the Sun and the donor star, respec- 
tively. 



2.3 Spin evolution of the NS 

In stellar evolution code, we also consider the spin evo- 
lution of pulsars as follows. With the spin-up of the NS, 
the accreting material would interact with the magneto- 
sphere of the NS. We simply define the magnetosphere ra- 
dius as the position that the ram pressure of the infalling 
material is balanc ed by the magnetic pressure of the NS 



Lam b et alJll973D. Und er assumption of spherical accretion 



Ghosh fc Lamblll979al lbl). the magnetosphere radius is 



= 1.6 x 10 s 



B s 



10 12 G 



4/7 



iMd 



lO^gs" 1 



-2/7 



(5) 



where B s is the surface magnetic field of the NS. Some obser- 
vations and analysis argued that the mass a ccretion of th e 
NS can lead to its magnetic field decay (s ee IWiiersHl997fl. 
Here w e adopt an empirical model given bv lShibazaki et al.l 
IHii ), i. e. 



Si 



1 + AM«/m B 



(6) 



where Bi is the initial magnetic field of the NS, AM aC c is 
the accreted mass of the NS, and me is ~ 1O _4 M0. 

When the NS rotation is too fast, the gravitational force 
of the accreting material at r m is less than its centrifugal 
force. The centrifugal barrier would eject the accreting ma- 
terial , and exerting a prop eller spin-down torque on the NS 
(111 arinov fc Sunvaevlll975l ). Namely, if the magnetosphere 
radius is greater than the co-rotation radius 



1.5 x 10 ! 



(Mns\ 
\M Q J 



1/3 



,2/3 



(7) 



where P B is the spin-period of the NS in units of second, the 
propeller effect occurs. The spin angular momentum loss 
rate via the propeller effect can be written as 



j p = 2MC[fi K (r, 



■0], 



(8) 



where S7K(V m ) is the Keplerian angular velocity at r m . When 
r m < r co , the accreting material is bound in the magnetic 
field lines to co-rotate with the NS, and is accreted onto its 
surface. Assuming rigid body rotation and the momentum 
of inertia / = 10 45 gcm 2 , the spin- up torque of the accreting 
material exerting on the NS is given by 



Jac = MacVGMNS-R, 



(9) 



where G is the gravitational constant, R is the radius of the 
NS. 

In addition, if r m is greater than the light cylinder ra- 
dius 



_ c _cP s 

ric ~ n ~ 2tt 



(10) 



the NS appears as a radio pulsar. As a result of magnetic 
dipole radiation, the spin angular momentum loss rate is 



3 RESULTS 

Based on the stellar evolution code and input physics de- 
scribed in Section 2, we calculated the evolution of large 
numbers of LMXBs. We take the initial mass of the donor 
star to be in the range 1.0 — 2.0 Mq (with a solar chem- 
ical composition Y = 0.28, Z = 0.02) 0, the initial mass 
of the NS to be 1.4 Mq. In addition, we assume that the 
NS is a dead pulsar, which already evolved to cross the so- 
called death line, and cannot radiate radio pulses. The spin- 
period of a dead pulsar s hould satisfy P 2 > B/1.7 x 10 11 
l|Bhattacharva et al.lll992T ). and the upp er-limit of the spin - 
period for normal radio pulsars is 11 s (jManchesterl 120041 ) . 
IWang et al. I (|201ll ) proposed that the minimum spin period 
of the NS is insensitive to its initial spin-period and magnetic 
field. Therefore, we adopt an initial spin-period of P s ,i = 10 
s and initial magnetic field of B = 10 12 G for the accreting 
NS. If the donor star evolves into a He white dwarf and the 
Roche lobe overflow ends, we stop the calculation. 

In Figure 1, we show an evolutionary example of an 
LMXB with an initial donor star of Ma,i = 1.5 Mq, and an 
initial orbital period of P or b,i ~ 1.5 day. Due to the loss of 
orbital angular momentum by magnetic braking, the orbital 
period decreases from 1.5 d to ~ 1.0 d before the mass ex- 
change. When the age is 2.15 x 10 9 years, the donor starts 
to overflow its Roche lobe, and mass transfer onto the NS 
commences. Because the material transfers from the more 
massive donor star to the less massive NS, the mass trans- 
fer firstly occurs on a thermal timescale at a high rate of 
~ 10~ 8 Meyr" 1 . At the same time, the surface magnetic 
field of the NS sharp decay to 10 9 G. In the initial mass 
transfer phase, the LMXB appear as a short-lived persistent 
X-ray source. With the decrease of the donor star mass, the 
mass transfer subsequently occurs on a nuclear timescale at 
a lower rate of ~ 10 -10 — 10 -9 Mq yr _1 . This rate is always 
less than the critical mass transfer rate 0- Therefore, the 
LMXB should be a transient X-ray source about 95% of all 
its life. After 1.24 x 10 9 years mass transfer, the NS grows to 
1.63 Mq. By gaining the material and angular momentum 
from the donor star, the spin-period P s of the NS continu- 
ously decrease to 1.8 ms. When the hydrogen-rich envelope 
of the donor star is exhausted, the degenerate He-rich core 
remains behind. The endpoint of the evolution is a BMBP 
consisting a recycled pulsar and a He WD with a mass of 
0.27 Mq, and with an orbital period of 27.3 d. 

The main aim of this work is to explore the progenitor 



4 Based on the detailed numerical c alculations for th e non - 
conservative evolution of close binaries. iTauris fc Savoniid jl999t) 
concluded that binaries containing a donor star with mass of 
> 2.0 Mq would experience a common envelope evolution, and 
evolve into BMSPs with a short orbital period (< 10 days) and a 
heavy CO or ONeMg white dwarf. 

5 There exist a lot of spikes and dips in the mass transfer rate 
curve, which are the results that the donor star exhau sted its core 
hydrogen and deviated thermal equilibrium JlH |2004) . 
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Figure 1. Evolutionary track of an LMXB with ; = 1.5 Mq, and P or b i = 1.5 day, which can evolve into a BMSP. The solid and 
dotted curves represent the evolution the NS mass and the magnetic field in the left panel, the donor star mass and the orbital period 
in the middle panel, and the mass transfer rate and the spin period of the NS in the right panel, respectively. 
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Figure 2. Distribution of the initial orbital periods P rb,i and 
the initial donor star masses M^j of LMXBs that can evolve into 
BMSPs via the recycling evolutionary channel. 



properties of BMSPs formed by the recycling evolutionary 
channel, therefore we have calculated the evolution of large 
numbers of LMXBs with different initial orbital periods and 
donor star masses. In Figure 2 we present the progenitor 
distribution of BMSPs in Ma,i — Porb.i diagram. The regions 
enclosed by the solid, dashed, and dotted curves represent 
the distribution areas of LMXBs that can result in a BMSP 
with a spin period of 10 ms, 20 ms, and 30 ms, respectively. 
Our results show that all NSs in LMXBs have a chance to be 
spun up to millisecond period, and the final fate strong de- 
pend on the separation of the binary. When the initial mass 
of the donor star is located in the range of 1.3 Mq — 1.6 Mq, 
the initial orbital period have a wider distribution from 1.0 
day to 20 days. When the initial mass of the companion is 
between 1.0 Mq and 1.4 Mq, the system cannot produce a 
BMSP unless the initial orbital period is less than 2.0 days. 



Figure 3. Predicted relation between the orbital period P or b 
and the white dwarf mass M^d for low-mass binary pulsars. 
The filled circles, the solid curve and the d otted curve denote 
our ca lculat ed results, the relation o btained bv lTauris fc Savoniiel 
lll999h . and lRappaport et al.l lll995h , respectively. 



Beyond these areas, BMSPs cannot be formed due to ei- 
ther a low spin-up efficiency or unstable mass transfer. For 
donor stars with a mass of 1.4 -1.6 Mq, a lower mass accu- 
mulation and spin-up efficiency of the NS result in an upper 
limit of initial orbital period. However, for massive donor 
stars with a mass of 1.7 -2.0 Mq, the upper limit on the 
orbital period originates from the dy namical instability of 
mass transfer (|Willems fc Kolbl 120021. In particular, m our 
calculated grids there exist two binaries which can produce 
a sub-millisecond pulsar with a spin period of 0.9 ms. Both 
LMXBs have a donor star with an initial mass of 1.9 - 2.0 
Mq, and an orbital period of 1.2 — 1.3 days. 

The evolutionary endpoint of most LMXBs is BM- 
SPs consisting of a millisecond pulsar and a low-mass He 
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Table 1. Observed parameters for 17 low-mass binary puslars. 
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Figure 4. Distribution of BMSPs in the spin-period P s of BM- 
SPs vs. the final orbital periods P or b diagram. The open circles, 
solid circles denote our calculated results and the observed data, 
respectively. 



white dwarf, which is called low-mass binary pulsar (LMBP) 
l|Stairs!l2004 iTauris fc van den Heuvelll2006l ). Stellar evolu- 
tion theory predicts a tig ht relation betwe en the core mass 
of giants and their radius (|Joss et al.lll987l ). During the evo- 
lution of LMXBs, the giant should overflow its Roche lobe, 
and its radius relates to the orbital separation. When the 
giant envelope is exhausted, its core evolve into a white 
dwarf. Therefore, the final orbital period of LMBPs should 
be correlated with the mass of the wh ite dwarf companion 
(|Savoniie f l987; Ra ppaport et al.l ll995'). Previous works pre- 
sented a simple relation between the orbital period P or b and 
the white dwarf mass Mwd for l ow-mass binary pul s ars (s ee 
also iRappaport et al. I { 19951 ) and lTauris fc Savoniid |l999j)). 
In Figure 3, we show our obtained low-mass binary pulsars 
by filled circles in P or b — A/wd diagram. It is clear that our 
ca lculated results are consi stent with the relation obtained 
bv lTauris fc Savoniid (| 19991 ). To compare with observations, 
we summarize the observed parameters for 17 low-mass bi- 
nary pulsars in Table 1. In Figure 4, we compare the calcu- 
lated results with the observed data in the P or b — Ps plane. 
It seems that our evolutionary model can account for the 
formation of part BMSPs. However, it is difficult for our 
evolutionary scenario to produce BMSPs with a short spin- 
period (3-8 ms) and a long orbital period (> 70 — 80 day). 

In Figure 5, we show the distribution of the final ac- 
creted mass and the final spin-period of NSs. One can see 
that, if NSs accretes a mass of > 0.1 Mq, they can be spun 
up to < 10 ms. In our calculated results, there exist 3 NSs 
that can accrete mass of > 0.6 Mq. Recent Shapiro de- 
lay measurements of PSR J1614-2230 suggested that it is 



a massive MSPs (~ 2 Mq), and with a CO white dwarf of 
~ 0.5 Mq l|Demorest et alfeoiOT lFl. We expect the discovery 
of LMBP with a massive NS like PSR J1614-2230 to test our 
evolutionary results. 



6 Recently, iLin et all l l201ll) proposed that this MSP may origi- 
nate from an IMXB with massive NS of 1.6 M q. Another work 
performed bv lTauris. Langer fc KrameJ l l201ll ) also support the 
viewpoint that the NS in PSR J1614-2230 was born massive. 
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Figure 5. Distribution of our simulated results in the spin-period 
P a of BMSPs vs. the accreted mass of NSs AM^jg diagram. 



4 DISCUSSION AND SUMMARY 

Using a stellar evolution code, in this work we have inves- 
tigated the formation of BMSPs formed by the canonical 
recycling evolutionary channel. In calculation, we take into 
account the influence of thermal and viscous instability of 
an accretion disk and propeller effect on the mass transfer 
process and the spin evolution of NSs. Orbital angular mo- 
mentum loss mechanism by magnetic braking, which orig- 
inate from the coupling between the stellar winds and the 
magnetic field of the donor star, is also considered. Our main 
results are summarized as follows. 

1. We obtain the initial donor star mass and orbital 
period distribution of the progenitors of BMSPs formed by 
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the recycled channel (see Figure 2). Our results show that 
all LMXBs with a donor star of 1.0 - 2.0 Mq have a chance 
to evolve into a BMSP with a spin-period of P s < 10 ms. 

2. The final fate of LMXBs strongly depend on the ini- 
tial donor star mass and the separation of the binary. When 
the donor star mass is in the range 1.0 - 1.4 Mq, only LMXBs 
with a short orbital period (P or b < 2.0 day) can evolve into 
low-mass binary pulsars. However, for a higher mass donor 
star, the upper limitation of the orbital period that can re- 
sult in birth of binary millisecond pulsars is 2.0 - 3.6 days. 

3. Our calculated results show that, if the NS accretes a 
mass of > 0.1 Mq, it can be spun up to millisecond period. 
In addition, it is possible that few MSPs gain a mass of 

> 0.5 Mq. 

4. It is difficult for our evolutionary scenario to produce 
a sub-millisecond pulsa rs. This result is consistent with th e 
conclusion obtained bv lFerrario fe Wickramasinghd l|2007l ). 
However, AIC evolutionary chan nel may produ ce a sub- 
millisecond pulsar (or quark star) l|Du et al.ll2009h . 

Obviously, our evolutionary results depend on the pa- 
rameterized input physics, especially the magnetic braking 
model, the duty cycle, and the magnetic field decay model, 
which have not been fully understood. Firstly, the loss of or- 
bital angular momentum plays a vital role in the evolution 
of LMXBs, hence magnetic braking model can influence the 
final orbital period of BMSPs. Secondly, the duty cycle can 
influence the outburst timescale and the mass growth of the 
NS. Therefore, a large duty cycle can result in the birth of 
MSPs with a short spin-period. In addition, the duty cycle 
may relate to system parameters (Lasota 2001), and may 
also evolve with the orbital period and mass transfer rate. 
Thirdly, in our input physics the magnetosphere radius re- 
lates to the field decay model, while their relation is not 
sensitive. Some uncertainties mentioned above may be re- 
sponsible for the discrepancy between our simulated results 
and observational data in Figure 4. A large duty cycle and 
a weak magnetic braking model may produce BMSPs with 
a short spin period (3-8 ms) and a moderate long orbital 
period (> 70 — 80d). Certainly, if our evolutionary model is 
correct, there may be other evolutionary channel to BMSPs 
such as AIC process of massive white dwarfs. 
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